Transverse momentum dependent distribution functions 
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Transverse parton momentum dependent distribution functions (TMDs) of the nucleon are studied in 
a covariant model, which describes the intrinsic motion of partons in terms of a covariant momentum 
distribution. The consistency of the approach is demonstrated, and model relations among TMDs are 
studied. As a byproduct it is shown how the approach allows to formulate the non-relativistic limit. 
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I. INTRODUCTION 

Studies of hard scattering processes such as inclusive deep-inelastic lepton nucleon scattering (DIS) have given 
rise to a good understanding of parton distribution functions, which tell us how the parton momenta parallel to 
the nucleon momentum are distributed. One way to gain insights into the partonic quark-gluon substructure of the 
nucleon beyond this one-dimensional picture is to consider transverse momentum dependent ('unintegrated') parton 
distributions (TMDs) %. These objects can be accessed by observing transverse momenta of, e.g., hadrons produced 



in semi-inclusive DIS 
factorization 
first data 
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SIDIS) or dileptons produced in the Drell-Yan process 
13l Il4fl . Much theoretical progress was made lla, [la LL7J, LL 



HT5gj54TT5^ 





thanks to 
MEL and 



give rise to phenomenological insights 



Nevertheless, presently model studies [gl SlM tsl HUM S^SSIzlSEIzlIzi eI^Y an important 
role. It is worth to recall that important insights concerning the very existence 9j or universality [161 ] of effects were 
made on the basis of model studies, see [79J for a review. Moreover, model results help to sharpen our physical 
intuition on these novel objects, and can be used to make estimates for planned experiments. Another aspect is that, 
thanks to the far simpler dynamics as compared to QCD, one may find relations among the different TMDs in some 
[M [M |M [M [6?J though not all [M [Ml models. As all TMDs are a priori independent structures, any such model 
relations among TMDs are not expected to hold in QCD. 

It is interesting to ask, however, whether such relations could neverless be satisfied in nature at least approximately. 
Recalling that the nucleon is characterized by 8 leading-twist 0] and 16 subleading-twist [23[ TMDs, such approximate 
relations could be valuable, for the interpretation of first data, or for estimates for new experiments (25j . 

In order to judge to which extent a particular model relation among TMDs might be respected by nature, till 
we know the answer from experiment, it is helpful to understand under which general conditions in a model this 
relation holds. For example, suppose a model relation relies on the SU(6) spin-flavour symmetry of the nucleon wave- 
function. We know from experiment that the SU(6) symmetry concept is useful — within certain limitations [78ll80l|. 
This implies that model relations based on SU(6) symmetry are respected in nature, if at all, at best within similar 
limitations. It happens that all model relations among TMDs observed so far have been found in models based on 
SU(6) symmetry [M, H| . 

The purpose of this work is to study TMDs in the covariant model of the nucleon proposed in Ref. [8l[ which makes 
no use of SU(6) spin- flavour symmetry. Some of the results presented here were discussed previously in [66| . In this 
model the intrinsic motion of partons inside the nucleon is described in terms of a covariant momentum distribution. 
The model was applied to the study of unpolarized and polarized parton distribution functions accessible in DIS 
fi(x), gf(x) and g^(x), and extended to c omp ute the transversity distribution hf(x) [82l. l83l l84l. l85l |86| . 

In this work we generalize the approach |8ll , l82l [H, [U, [85], [8(| to the description of TMDs. In particular, we focus 
on the so called T-even, leading twist TMDs, and pay particular attention to the demostration of the consistency of 
the approach. We shall see that certain model relations among TMDs hold even without invoking SU(6) symmetry. 

This note is organized as follows. In Sec. [IT] we briefly introduce TMDs. In Sec. IIIII we introduce the model, 
and review previous works. In Sec. [TV] we generalize the approach to the description of TMDs, and in Sec. [V] we 
demonstrate its consistency. In Sec. I VII we discuss the model relations among the polarized T-even TMDs. In Sec. I VIII 
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we apply the approach to a study of TMDs in the non-relativistic limit, before we summarize and conclude in Sec. lVIlH 
The Appendices contain details of the calculations, and supplementory results. 



II. TMDS 



In this Section we introduce and define briefly TMDs. With the use of light-cone coordinates, a ± = (a ± a 1 )/-^, 
TMDs arc defined in terms of light-front correlators as 



dz d 



(2tt) 3 



,i P z (7V(Pj (Q) W(Q) Z) path) ^ ( z ) |iv(F, S)) 



(1) 
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In SIDIS the singled-out space-direction is along the momentum of the hard virtual photon g M = ft? , q 1 , 0, 0), and 
transverse vectors like pr are perpendicular to it. The path of the Wilson- link depends on the process 
the nucleon rest frame the polarization vector is S = (0, — Sl, St) with S% + S T = 1. The negative sign in front of 
Sl is because by convention [24| the nucleon has positive helicity, i.e. Sl > 0, if it moves towards the virtual photon. 

The information content of the correlator (fTJ) is summarized by eight leading-twist TMDs Q, that can be projected 
out from the correlator ([!]) as follows (for convenience we will often suppress flavour indices) 



f 



tr 



7+ (f>{x,p T ) 



tr 



7+75 0(x,p T ) 



tr 



\x,Pt, 



h 

= S L gi 



Sip h\ 



s^PtSt , x 
Mn hT 

Pt ■ St ± 
■9it 



S L 



N 

Pt 



{PtPt-iPt^St 



N 



L 1L 



M 2 

N 



e jk p T 



M 



N 



(2) 
(3) 
(4) 



where the space-indices j, k refer to the plane transverse with respect to the light-cone, and e 32 = — e 23 = 1 and 
zero else (which is consistent with e 0123 = 1). Integrating out transverse momenta in the correlator (JT|) leads to the 
three 'usual' parton distributions known from collinear kinematics (x) = J d 2 p r jf(x. p^) with j = f, g, h f87L fS8j] . 
Dirac-structures other than that in Eqs. ((21 [21 SI) lead to subleading-twist terms [23l l24j|. 



III. THE COVARIANT MODEL OF THE NUCLEON, AND ITS APPLICABILITY 

In this Section we first briefly introduce the model, and sketch the calculation of the 'collinear' parton distribution 
functions done so far, namely fi(x), gl(x), g T (x) and h^(x). Then we discuss the applicability of the approach to 
the calculation of TMDs which will be done in Sec. IIVI 

The starting point for the calculation of the chirally even functions accessible in DIS, fi(x), gi(x), g T {x) [HI, [H, [83[ , 
is the hadronic tensor. The latter is evaluated in the Bjorken-limit, i.e. in the limit that the four-momentum transfer 

from the lepton beam to the nucleon with momentum P is such that Q 2 = —q 2 and Pq — ► oo while x = Q 2 / (2Pq) 
is fixed. In the model it is assumed that unpolarized DIS can be described as the incoherent sum of the scattering 
of electrons off non-interacting quarks, whose momentum distributions inside the nucleon are given in terms of the 
scalar function G(pP/M). Here p and P are the momenta of the quark and nucleon, and M is the nucleon mass. 
Though all expressions can always be formulated in a manifestly covariant way, it is convenient to work in the nucleon 
rest-frame, where the momentum distribution becomes G(p°) with p° = \J p 2 + m 2 . Here m denotes the quark mass. 
Clearly, the distribution of the quark momenta in the nucleon rest frame is rotationally symmetric. 

Applying these ideas to the description of the symmetric part of the hadronic tensor has shown that in the model the 
Callan-Gross relation among the unpolarized structure functions holds exactly, and the unpolarized parton distribution 
function is given by [8l| (notice that G(p°) depends on flavour, which we suppress for brevity) 

ff(x) =J^ G(p°) S (l^f- - «) (pP - p 1 ) . (5) 

Next we review the calculation of g1(x) and g T (x) appearing in the anti-symmetric part of the hadronic tensor. 
For a single quark the latter would be given by W^ q = •mt a ^^ v q >1 w v , where q^ is the momentum transfer from 
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the electron in DIS, and w u denotes the polarization vector of the quark. In the model — assuming the covariant 
distribution of polarized quarks to be given by H(pP/M) — the anti-symmetric part of the hadronic tensor of the 
nucleon is given by 



W, 



^ H(pP/M)S((p + qf m 2 )W^ § ^~ q e a ^ 



M 



(6) 



The second expression in ([5]) is given in the nucleon rest frame choosing g M = (g°, g 1 , 0, 0) and holds in the Bjorken 
limit (more precisely: here and in the following for the steps marked by 'Bj' the condition Q 2 ^> 4M 2 x 2 is essential). 

Notice that the covariant distribution of polarized quarks can be expressed as H(p°) — G + (p a ) — G~(p°), where 
the indizes (±) refer to the respective quark polarizations which are parallel (+) or anti-parallel (— ) to the nucleon 
spin. In this notation the covariant distribution of unpolarized quarks in ( Op i s G(p°) — G + (p°) + G~(p°). 

The most general expression for the covariant quark polarization vector [82[ is given by 



pS 



pP + mM 
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where 5 M denotes the nucleon polarization vector given in the nucleon rest frame by S 1 ^ = (0, S) with l^l = 1. The 

tA 



evaluation of Eqs. (J6j [TJ) and comparison to the general Lorentz-decomposition of W^g, namely 



Pq 



9i 



(Pq)S v - {Sq)P» 
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.92 



yield the following results for g\{x) and g T (x) = gf(x) + g\ (x) 



9t( x ) 
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In the model the Burkhardt-Cottingham sum rule [90] is satisfied. When neglecting terms proportional to m also the 
Efremov-Leader-Teryaev sum rule [9l| holds, while g T (x) is given by the Wandzura-Wilczek (WW) approximation [92j 



g T (x) 



WW 



(ii) 



as was proven in [82 . [83j . Notice that in QCD what is neglected are not only mass terms but also pure twist-3 terms 
[92| . That in the model such pure twist-3 ('interaction dependent') terms in gT( x ) are absent, is consistent because 
in our approach the quarks are assumed to be free. 

The chirally odd transversity distribution function cannot be accessed through the hadronic tensor and DIS. For 
theoretical purposes, however, one may consider the auxiliary polarized process described by the interference of a 
vector and a scalar current. On the quark level this interference is described by — e a i3\vP^q*w 1 ' from which one 
obtains — in analogy to the procedure in Eq. ([6]) — the following expression for the nucleon 



d 3 p 



H(pP/M)S((p + q) 2 - m 2 )Tl 



2p- q e ^ q J "p0~ H{p )5 {-^r- - T 



Pp/M " yt " '"' '"^ ' ~ 1 ' "" '" a 2Pq 
The general Lorentz-decomposition in this case reads (J — 2, 3 is a 'transverse index' with respect to q and P) 

2M(-iy a T a = S T h\{x) 

one obtains after evaluating (fT2"|) with the following result for the transversity distribution function [8~ij 



h\{x) 



d3 P tt 1 0^r/P° 
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(14) 
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We remark that in Eqs. ([51 [TUl [T^|) we did not distinguish momentum distributions in differently polarized nucleons. 
In general one might suspect the covariant distributions to be different. In QCD, if nothing else, different evolution 
properties clearly distinguish chirally even vs. odd, and twist-2 vs. twist-3 distribution functions. In the model, 
however, it is natural to assume the distributions in longitudinally and transversely polarized nucleons to be equal. 
In order to understand that this assumption is indeed natural, we recall that the approach is covariant. Therefore 
one may go to the nucleon rest frame, where it certainly makes no difference whether the quarks in the nucleon are 
polarized longitudinally or transversely (with respect to the space component of the four- vector q in DIS or SIDIS). 1 
Since in this model the quarks are non-interacting, it does not matter how they are polarized — because, for example, 
there are also no spin-orbit- or spin-spin-interactions. (Sec. IVI A"l will show that the covariant distribution in various 
polarized TMDs must be equal, in order to comply with QCD.) 

However, this by no means implies that the parton distributions describing longitudinally and transversely polarized 
quarks, gi(x) and h\{x), are equal. They are, in fact, rather different even if described in terms of the same covariant 



momentum distribution H(p°) 8J|. By introducing adequate normalizations (as dictated, e.g., by SU(6) symmetry) 



one could furthermore relate the polarized covariant distribution function H(p°) to the unpolarized one G{p°). This 
is, however, a severe restriction and simplification of the model, which we do not need in general. 

When extending the approach below to the description of TMDs it is important to keep in mind the following point. 
The QCD definition of a parton distribution function includes a Wilson line, which in DIS describes the interaction 
of the struck nucleon with the target remnant. It is possible to find a gauge in which the Wilson line drops out, and 
the partons seem 'non-interacting' — an idea eventually underlying the parton model in general, and the approach of 
Refs. l82l. l83l [84l [85l [8(| in particular. When dealing with TMDs, however, the Wilson line cannot be 'gauged away' 

In the present framework we have no tool to include the effects of the Wilson line, and therefore the description of 
the so-called 'naively T-odd' parton distribution functions, the Sivers function and the Boer-Mulders function hj~, 
is beyond the scope of the approach. These TMDs crucially rely on the initial- and final-state-ineractions encoded in 
the Wilson line [1, [l(| , and are expected to be absent in our framework. 

Finally, we remark that the model is 'opposite' to the Gaussian ansatz for TMDs in the following sense. In the 
Gaussian ansatz one assumes the extreme situation that distribution of longitudinal momentum (i.e. ^-dependence) 
and the distribution of pt are decoupled. For example, one has fi(x,p^) = f\{x) exp(— p 2 -. / {p 2 ^)) / {-k {p\)) . It is even 
possible to "switch off" pT-effects: in the limit (p\) — > one has fi(x,p^) = fi(x) S^ 2 >(pr). In contrast to this in the 
present model the longitudinal and transverse motions are coupled "maximally" . It is not possible to "switch off" 
PT-effects. As a consequence one has, e.g., interesting implications for the quark orbital motion (86j . 



IV. EXTENSION OF THE APPROACH TO TMDS 



In this Section we extend the approach to the description of TMDs. Since none of the new TMDs in Eqs. (J2H4J) is 
accessible directly via the hadronic tensor in DIS or via the auxiliary process explored for the calculation of h\(x), we 
need to establish a more general relation in the model to the correlators ((2HU). For that we observe that the model 
expressions for the anti-symmetric part of the hadronic tensor © or the auxiliary current (|12[) include integration 
over d?p = dp 1 d 2 pr with d 2 pr = dp2dp%. In the following we will explore the consequences of what happens if one 
does not integrate out transverse momenta in these expressions, and demonstrate the consistency of this approach. 

The 'integrated' symmetric part of the hadronic tensor, to which ff is related, was studied in Ref. 81]. The study of 
its 'unintegrated' version would give model results for and the T-odd /{^?, as revealed by the correlator in Eq. ^fy. 
However, the description of the Sivers function is beyond the scope of our approach, see Sec. IIII1 and we therefore 
start with the more interesting case of the correlator @ which describes gf(x,pT) and the T-even TMD gf T (x,pT)- 
(We shall come back to ff at the end of the next Section.) 

In order to access the information contained in the correlator §3§ we consider the transverse space components 
(j, k = 2, 3) of the 'unintegrated' anti-symmetric part of the hadronic tensor in Eq. ©. We work in the nucleon rest 
frame with the nucleon polarization vector as introduced in the sequence of Eq. ((T|) and choose = (q°, q , 0, 0). 



Since is space-like q 2 < 0, its space-component q is non-zero in any frame, and always provides an axis for the quantization of the 
nucleon spin. 
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Then, using ([7]) we obtain in the Bjorken-limit 



-S L pV-p 1 



P' 
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PtS t P° -p 1 

2 p° + TO 



(15) 



We recognize two contributions in (|15[) . one proportional to the longitudinal nucleon spin component Sl and one 
proportional to the projection of the nucleon spin on the transverse parton momentum. These contributions coincide 
exactly with the decomposition of the correlator in Eq. ([3]). Thus, from the comparison of the coefficients we obtain 
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(16) 
(17) 



Notice that there is no arbitrariness concerning an overall prefactor, because in the integrated case we reproduce 
2MWjf.(x) = —€jk^L9i(x) in agreement with the general Lorentz-decomposition for the transverse components of 
the anti-symmetric part of the hadronic tensor. The other components of W^ v describe subleading twist structures, 
for example g^(x) in the integrated case, which we do not consider in this work. 

Now we wish to access the information content described in the chirally odd correlator (j4|). For that we consider 
the 'unintegrated' version of the auxiliary current T a in Eq. (|12[) . For T a (x,pr) contracted with e ja we obtain: 



2.u i -n.-'-'T,A.r.,r 1 i = / ^H(p )S fl ' [ 1,1 
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M J \ 1 ^ 1 p° + m ± p u to 

In order to easier compare to Q we rewrite the decomposition of that correlator as often done [f| 0] as follows 
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where we suppressed the arguments x, px of the TMDs for brevity. By comparing the coefficients in Eqs. ([TBI [T9"]) we 
read off the following results: 



h q 1T (x,p T ) = 
h il{x,Pr) = 
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hi q (x,pT) = 
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CONSISTENCY OF THE APPROACH 



It is necessary to demonstrate the consistency of our approach. For that we remark first that by integra ting the 
expression for g1(x,pT) in Eq. (|16p over transverse momenta we recover the model result for gf(x) derived in [82|and 
quoted in Eq. Next, by exploring the connection of transversity to the functions Iiit(x,Pt) and h 1 L T (x,pT) [61,0] 



hi(x,p T ) = h 1T (x,p T ) + h XT {x,p T ) 
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and integrating over transverse momenta we recover the correct result for transversity derived in [84j and quoted 
above in Eq. (fl4|) . This means that the Lorentz-decomposition of the structure (| 18|) in the model is consistent with 
the Lorentz-decomposition of the correlator in Eqs. (fH [TO]). Since the model is covariant, this is an expected feature. 
We finally notice that there is no polarization-independent term in (fl8|) meaning the absence of the Boer-Mulders 
function — as expected in the present framework, see Sec. IIIIl 

We observe that the results for g q (x,pr) and h\{x,pr) derived here could have been simply 'guessed' from the 
results for the integrated functions gf(x) and h\(x) by 'skipping' the integration over d 2 pr- This is by no means 
trivial, because in the respective TMDs there could have been structures — e.g., of the type F(p° ,p x )[p2 ~ P3) with 
some function F(p°,p 1 ) making the integrals converging — which would drop out in the expression for the integrated 
functions due to rotational symmetry in the transverse plane. On the other hand, in the present approach described 
by a free Hamiltonian (that commutes with the momentum operator) the model expressions for TMDs could be 
written as nucleon expectation values of certain polynomials of the momentum operator. Then it is clear that such 
'multipol-terms' in TMDs are forbidden by the Wigner-Eckart theorem. After these considerations we conclude that 
the model expression for the unintegrated unpolarized distribution function is given by 



ff(x, PT ) = / ^ G(p°) S 



'<■'"' - x ) {p v - p ) 



M 



(25) 



As a next important consistency check let us test whether the model results satisfy positivity constraints, and for 
that we need the expression for the unpolarized distribution function (|25p . since the inequalities for TMDs we wish 
to verify read [H| 



\h\(x,pr)\ < 



,-L(i)g 

'IT 
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Pt 
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fi(x,PT) 2 - g q (x,p T ) 2 
f?(x,p T ) 2 - g q {x,p T ) 2 



where the 'unintegrated' transverse moment of a TMD is defined as 



? -(D 
Ji 



(x,Pt) 



Pt 
2M 2 
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(26) 
(27) 
(28) 
(29) 

(30) 



Notice that the T-odd functions f^^ix^pr) and hf^ q \x,pt) are absent in our approach. A direct test of the 
inequalities is actually difficult because different covariant functions G(p°), H(p°) apprear in the TMDs. One way to 
proceed is to assume SU(6) spin-flavour symmetry of the nucleon wave function. If one assumes SU(6) symmetry, the 
inequalities (|26H29p are manifestly satisfied in our framework, see App. [X] If one does not, the positivity conditions 
(I26H29I) are 'translated' into certain constraints among the covariant momentum distributions G(p°) and H(p°), see 
Ref. \8A\ , where the pr-integrated version of (|26| known as Soffer bound [93j was discussed in this way. 

We conclude this Section with the observation that so far our approach satisfied all imposed consistency checks. 
Further consistency tests will be provided in the next Section, where we shall see that our approach satisfies certain 
exact relations as well as model relations among different TMDs found also in other relativistic quark models. 



VI. RELATIONS IN THE MODEL AMONG TMDS 



It has to be stressed that in QCD all TMDs are independent structures, and it is not possible to find exact relations 
among them that would allow to express one TMD in terms of other TMDs. However, especially in models without 
gauge field degrees of freedom [H,[64, 65], it might be possible to find relations among different (T-even) TMDs. Such 
model relations are of interest by themselves, and might be supported by data within certain 'model accuracies'. It 
would be interesting to know the general conditions a quark-model must satisfy in order to fullfill such relations. 
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In order to recognize more easily the relations among the different T-even TMDs let us introduce the following 
compact notation for the measures 



{dp} = — 6 



p° — p 1 



u _ dp 1 H(p°) f P ' 



WP > = ST W> 



p u p u + m 



M 



(31) 
(32) 



The measure (|5T|) is positive definite, while the sign of the measure (j3"2"|) depends on the sign of H(p°). Then the 
various TMDs can be written as follows 
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(33) 
(34) 
(35) 
(36) 
(37) 
(38) 



where we remind that p° = \J p\ + p£ + m 2 and p° — p 1 = xM due to the delta-function in the measures. When 
deriving the second equality in (|35[) one may make use of the identity p^ = xM(p + p 1 ) — m 2 valid due to the delta- 
function. From the expressions (|34H38|) we can read off numerous model relations among polarized TMDs, which we 
shall discuss in the following. 

Notice that none of the relations discussed in the following involves the unpolarized parton distribution function. 
Such relations are impossible in our approach, simply because there is in general no way to connect the different 
covariant distributions G(p ) and H(p®). If (and only if) one makes an additional model assumption, namely assumes 
the SU(6) spin-flavour symmetry, then we obtain relations including fi(x,pr) and other TMDs, see App. 151 



A. An exact relation in QCD, and its consistent realization in the model 



Let us first discuss an exact relation which is valid in the model and which involves g^(x). (Later we shall discuss 
also several approximate relations involving this twist-3 parton distribution function.) For that we introduce the 
measure {d 3 p} = {dp 1 }d 2 pT which allows us to rewrite the model expression (jlOp as 



# (s) = J {d 3 p} 



m(p° 



1 -2 

2 PT 



From Eqs. ([35l [36l [39]) we find that the following QCD relation [24[ is satisfied consequently in the model 

x g q T (x) = g^ 1)q {x) + jjK {x) 




(39) 



(40) 



This further demonstrates the consistency of our approach, although g^(x) ^ in general, because in our model such 
pure twist-3 terms (quark-gluon-correlations) are consequently absent. We learn two further important lessons. 

First, the relation (|40[) crucially relies on the fact that h\(x), g x ^ are described in terms of the same covariant 
distribution function H(p°). In other words, to be consistent with QCD, we actually have no choice but must work 
with the same covariant distribution function H(p ) for all polarized functions (c.f. the discussion in Sec. lIIII) . Second, 
we clearly see that the model parameter m is really to be identified with the current quark mass in QCD. 
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Let us remark that in QCD the relation (|4T))) is valid also in 'unintegrated version', i.e. with the TMDs not integrated 
over px- Here we confine ourselves to the 'integrated relation' (|40[) . as we have not derived the model expression for 
g^{x,pT) (though, in view of the experience with g\ and h\, presumably it is given by (fTUl with p^-integration 
omitted and the 'unintegrated version' of ([40)1 is valid, too). 



B. Exact relations among leading- twist TMDs in the model 



Next we focus on a class of relations among leading twist TMDs which are exact in our approach, and can hold in 
models only. By comparing the expressions in Eqs. ([361 I3T[) we observe the following exact relation in our model 



g^(x,p T ) = -h^l(x,p T ) ■ 

This relation was observed previously in the spectator model of Ref. [63[ and the constituent model [6£ 
Next, by comparing Eqs. ([Ml 1551 155)) we see that the model results satisfy the relation 



gf(x,p T ) - h\(x,pr) 



^(1)9 
IT 



{x,p T ) 



(41) 



(42) 



This relation was first observed in the bag model [6J| and it is also valid in the spectator model of Ref. [63|]. It 
was argued [11] that (|42|) could be valid in a larger class of relativistic models. It is thus gratifying to observe 
that subsequently (|42p was confirmed in the relativistic constituent quark model |65| ] and now also in our approach. 
Interestingly, the relation (|4"2"]) is not supported in the spectator model version of Ref. [68| . 

Both quark model relations, Eqs. (|4T1 [42]) . are not supported in models with gauge-field dregrees of freedom [691 ] . 
This observation is in line with the expectation that even if the relations ([4T1 l42jl were valid in QCD at some scale 
(which, of course, does not need to be the case) they would be spoiled at any different scale by evolution effects that 
clearly discriminate chirally even and odd functions. 

Finally, from Eqs. (|20l [2X1 [22)) we find the following remarkable exact relation 



h\(x,pr) hf-£{x,p T ) , 



(43) 



that was not observed before in literature to best of our knowledge and connects only chirally odd TMDs — in contrast 
to (|41[ |4"2"[) . This non-linear relation is not obeyed in the spectator model [63[ . Combining (|4"T1 |4"3"[) we find 



9it( x >Pt) 



h\{x,p T ) h^(x,p T ) 



(44) 



which again mixes chirally odd and even TMDs (though the product of two chirally odd objects 'conserves' chirality). 
From Eqs. ([3X1 [3"5)) we find also relations among the signs of TMDs, for example: 



signal) = signfo?) , 
sign(/i?) = - sign(Zi^) 



Notice that ([46)1 could be concluded as a corollary from the result (|43[) . Also from Eqs. 
([43} and 03}, we find 

\h\{x,p T )\ > \g\{x,p T )\ for p T > , 



(45) 
(46) 

or from combining 
(47) 



i.e. the modulus of the transversity distribution function is larger than that of the helicity distribution function. This 
inequality survives integration over px and x, and we obtain for g^ = J &xh\(x) and g\ = J dx g\{x), the tensor and 
axial charges, the relation 



\9t\ > \9a\ , 



(48) 



which was also observed in [84( , in many other models [9J, [95|, |96j , and in lattice QCD [97| • 

Although some of these results have been obtained before in literature, it is remarkable that in our approach all 
these relations follow without assuming SU(6) spin-flavour symmetry. 
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Relations in the model in the chiral limit 



From Eqs. (|35II38[) one can find further relations by observing that /if, g^, h^£, h-,^, are functions of the type 
const x (xM + m) n x /{dp 1 }, which looks 'so trivial' only due to the compact notation introduced in Eq. (|32|) . Recall 
that m is to be identified with the QCD current quark mass, see Sec. I VI A] whose effects are expected to be negligible 
in deeply inelastic reactions. We can formulate those relations in the chiral limit, and obtain 



(49) 

2h\(x,p T ) + xhH(x, PT ) =0(j^ , (50) 
hfl(x,p T )-xh±2(x,p T ) = o(j^\ . (51) 

Upon the use of (HI]) one obtains relations similar to ([5171 IB"!]) but with replaced by (— g^)- 

Remarkably, in the model it is possible to relate the transverse moments of the chirally odd TMDs hf£ and hfj, to 
the chirally even twist-3 parton distribution function g^(x) as follows 

h^ )q (x)+xg^x) = o(^\ , (52) 

hf^(x)+g T (x)=o(^\ . (53) 

Notice that (|52|) follows also from combining (|4"0"]) and (|4"Tj) . 



D. Wandzura-Wilczek (type) approximations 

We already mentioned the Wandzura-Wilczek approximation 92] which allows to connect the twist-3 distribution 
function g^(x) and the twist- 2 distribution function g\{x), see Eq. (|lip . In the model this approximation requires the 
neglect of quark mass terms. In QCD one has to neglect in addition pure twist-3 terms. 

An important practical application is that Eqs. (jlll I53[) allow to express the unknown h^}f (x) in terms of the 
well-known g\ (x) . In [661 ] we made use of this relation in order to estimate the transverse moment of pretzelosity, and 
used the results for predictions of SSAs in SIDIS. 

In a similar way, i.e. upon the neglect of pure twist-3 and quark mass terms, one obtains further 'Wandzura-Wilczek- 
type approximations' [25|, namely 

gtP q (x) « x t^-gl{y), (54) 
j x y 

htP*(x) « - x *f^-hl{y). (55) 

j x u 

Both Wandzura-Wilczek-type approximations are valid in our approach upon the neglect of quark-mass terms. The 
validity of (|5"4"]) follows directly from Eqs. (fTTl |4"0")) . The proof of the Wandzura-Wilczek-type approximation ({55]) is 
given in App. [C] 



E. Transverse momentum dependence 



The probably most exciting thing about TMDs is, of course, their ^-dependence. The power of the covariant 
approach with rotationally symmetric momentum distributions of quarks in the nucleon rest frame is based on the 
fact that this symmetry ultimately connects the distributions of transverse and longitudinal momenta. Some exciting 
consequences of this symmetry in the context of the spin content were discussed in [86j , and a detailed study of the 
effects of this symmetry in the context of TMDs is in preparation. 
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However, a couple of simple but already interesting conclusions on the 'mean transverse momenta' of TMDs can be 
drawn without modelling the covariant momentum distribution H(p°). Notice that all py-integrals given below are 
well-defined in our approach. 

Let us introduce the notion of mean transverse momenta moments of a TMD j\ as follows 

J dxj d z p T Ji{x,pt) 

and the n-th moment p-moment (p = \p\) of the covariant distribution function is defined as follows 

«p«)) = J d 3 pp n H(p°) . (57) 
Then we obtain the following relations valid in the chiral limit 

> = = (58) 

Bm (p|, 9 iS) = Hm(p|, hf D = \ ^ (59) 

J^>4*> = (60) 

It is instructive to learn that, although g\(x) and h\(x) are very different in the model 84], their mean transverse 
momenta coincide in the chiral limit. 



VII. NON-RELATIVISTIC LIMIT 



The assumption of a non-relativistic dynamics of light quarks in the nucleon is not realistic. Nevertheless certain 
conclusions from this limit are very popular — like, for example, the relation h\(x) = gi(x) which has often been used 
in literature to obtain order of magnitude estimates for effects of transversity. It is therefore worth to study how this 
limit can be formulated in our framework. This will yield non-relativistic limit results for TMDs. 

In the strict non-relativistic limit, we have particle conservation and the nucleon consists of exactly 3 (valence) 
quarks. Then we deal with the dynamics of constituent ('valence') quarks, whose momenta become negligible with 
respect to m, and whose binding energy becomes negligible with respect to the nucleon mass such that M = 3m up to 
relativistic corrections. The heavy constituent quarks obey spin-flavour symmetry which is introduced in the context 
of TMDs in App. [£]in Eq. ([51]) . 

The non-relativistic limit makes the following predictions for the collinear parton distribution functions 

Jim f«(x)=N q s(x-^) , limjl{x)= UmjUx)=P q 8(x-^\ . (61) 

If the momenta of quarks are not negligible with respect to their mass m, the 5-functions in Eq. (|61[) are spread out. 
The normalizations N q and P q in Eq. (|6ip dictated by SU(6) spin-flavour symmetry are given in Eq. (|Alj) in App. [A] 
One can check that all sum rules are correctly reproduced: 



dx fx(x) — N q ('normalization') (62) 
/ dxxf1(x) = 1 (momentum sum rule) 



dx (g?(x) - gt(x)\ = J dx (h?(x) - ht(x)\ = gf = gf 



5 

— (Bjorken sum rule/isovector tensor charge). 
3 
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The model is formulated in a covariant way, which means that it also can be applied to the situation when the motion 
of quarks is assumed to be non-relativistic, i.e. when \p\ <C m and po = m{l + 0(p 2 /m 2 )}. In order to formulate 
the non-relativistic limit in our model, we assume SU(6) symmetry and set the covariant momentum distributions in 
unpolarized or polarized nucleons equal, i.e. we assume G(p°) — > N q J(p) and H(p°) — > P q J(p) with J d 3 p J(p) = 1. Of 
course, J{p) strictly speaking depends only on the modulus \p\ (or p° = ^Jm? + \p\ 2 ) due to the rotational symmetry 
in the nucleon rest frame, but this notation is more convenient for the following. 

In the non-relativistic limit one expects only small momenta p — > to be relevant in the integral over d 3 p. Thus, 
it is natural to assume 



lim J(p) -» 8 {3) (p) 



(63) 



This is in some sense an 'axiom', and we have to verify that it yields to consistent results. For that we insert (f6"3"| in 
Eq. ([5]), and obtain (notice that p° — ► m for \p\ <C m) 



lim flix) = N q xM I -f S^(p) S 



d3 P rf3U^ s rfP°+P 1 



M 



xM 



s(™ 



N q 5 I x 



(64) 



Thus, our prescription (|63|) gives the correct non-relativistic result. What do we obtain for the collinear polarized 
distribution functions? Let us insert ([63]) in Eqs. ([9j [101 IT4]) We obtain 



lim g\{x) = lim g^(x) = lim h\{x) = P q 8 I x — — 



1 



(65) 



Thus, the polarized distributions <?i(x), g^(x), h\(x) become equal and correctly reproduce the non-relativistic result 
(|61[) . Thus, we see that our formulation of the non-relativistic limit also yields correct results for the polarized parton 
distribution functions. 

Let us now apply the non-relativistic limit to the description of TMDs. We obtain 



lim ff(x,p T ) = N q S[x- — 



\im^(x,p T ) = P q 5\^x 
lim h\{x,pr) = P q 8[x 

non-rcl V 

\\mg^(x,p T )= N c P q s(x 
lim hfji(x,pT) - 



N c P q 6[x 



N 2 



1 

N c 
1 

iV^ 
1 

No 
1 

No 
1 

iV^ 
1 



S (2) (pt) 
5 {2 \pt) 
5 {2 \pt) 
8 (2) {pt) 
5 {2 \pt) 
S (2 Hpt) 



(66) 
(67) 
(68) 
(69) 
(70) 
(71) 



In g^jl(x,pT) and H^(x,pt) the factors ^ = N c appear, while in h^(x,pr) the factor = ^f- appears with 
N c = 3 colours. These factors appear here somehow artificially because the nucleon mass was chosen in the Lorentz- 
decomposition of the correlators ([2H1]) to compensate the dimension of transverse momentum. Nevertheless, once 
one introduces the nucleon mass in this context (and sets the according 'units to measure' TMDs), the integrated 
functions g^(x), /i 1 ^(x), h^{x) are larger then the parton distributions g\{x) and h\(x). It even happens that 
\h^{x)\ > ff(x) as also observed in other models (641 . J6||. This is not in contradiction with positivity which 
constrains only the transverse moments of TMDs, see Eqs. (|27l l29|) . 

From Eqs. (|66H71[) we see that the transverse moments of all TMDs vanish. In particular, hfj <(x) vanishes. This 
is consistent from the point of view of the relation between helicity, transversity and (the transverse moment of) 
since in the non-relativistic limit, helicity and transversity distributions become equal (|65p . 



Ml 



pretzelosity, Eq. 

Thus, a non-zero transverse moment of pretzelosity [64j or any other TMD (as we learn here) can be considered to 
be a 'measure of relativistic' effects in the nucleon. Clearly, any effect of TMDs would dissapear from a cross section 
(or spin asymmetry). However, the TMDs themselves are all non-zero in the non-relativistic limit, see Eqs. (|66H7ip . 
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VIII. CONCLUSIONS 

We have generalized the covariant model developed in Refs. [HI, HH HH HH, HH to the description of T-even 
leading-twist TMDs in the nucleon. We have payed particular attention to the demonstration of the consistency of 
the extended approach. For example, we have shown that it g ives the familiar results for the 'integrated' functions 
known from studies of collinear parton distributions [8ll |82| . l83l . |84| , l85l . [86[ , proven that it satisfies inequalities among 
TMDs, and discussed that it yields results consistent with the large- N c limit, lattice QCD, and many other models. 

In particular, we have also shown that in the approach a relation, which is derived from the QCD equations of 
motion and connects several TMDs and a pure twist-3 ('tilde') function, is consequently satisfied in the model. In 
our covariant approach with free partons 'consequently' means that the 'tilde'-function is absent. 

We have rederived several known quark model relations among polarized leading- twist TMDs [H, [64], [6|| , and found 
several new relations so far not observed in models, without assuming SU(6) spin-flavour symmetry. In our approach 
these relations refer to a scale of several GeV 2 . Whenever previously such relations were observed, the corresponding; 
model explicitly made use of the SU(6) symmetry and the results referred to low hadronic scales [63|, HH, [65[. Not 
all quark models support these relations [68] , and by including gauge-field degrees of freedom [6^ such relations are 
definitely spoiled which one expects to be the case also in QCD. However, it remains to be seen whether in nature 
some of these relations could at least be approximately satisfied. 

We have also shown that the Wandzura-Wilczek-type approximation, which allow to approximate the transverse 
moments of gf£ and hf£ in terms of respectively g\(x) and h\(x) are valid in the model upon the neglect of quark 
mass terms. In QCD these relations are valid if one in addition neglects also pure twist-3 terms [25[. 

As an interesting digression, we have discussed how the covariant model framework can be used to formulate the 
non-relativistic limit for TMDs, and derived the non-relativistic limit results for all leading-twist T-even TMDs. In 
the non-relativistic limit all these TMDs are non-zero, however, their transverse moments vanish. Interestingly the 
non-relativistic approach is consistent with the basic features of the relativistic model calculations. 

In this work we focussed on the general aspects of TMDs in the model. Further consequences for TMDs due to the 
parton intrinsic motion, as well as phenomenological applications (see [66| for first results) will be discussed elsewhere. 

Acknowledgements. A. E. and O. T. are supported by the Grants RFBR 09-02-01149 and 07-02-91557, RF MSE 
RNP.2.2.2.2.6546 (MIREA) and by the Heisenberg-Landau and (also P.Z.) Votruba-Blokhitsev Programs of JINR. 
P. Z. is supported by the project AV0Z10100502 of the Academy of Sciences of the Czech Republic. 

APPENDIX A: PROOF OF INEQUALITIES 

In this Appendix we prove that the inequalities ([27H29]) are satisfied, if one assumes SU(6) symmetry and if one 
assumes all TMDs to be described in terms of the same covariant distribution J(p ) normalized as J d 3 pJ(p°) = 1. 
Then, in SU(6) with N c = 3 denoting the number of colours [98(, the TMDs of definite flavour are given by 

gl(x) = P q gi(x), P u 

The 'flavour-less' functions introduced in (|A1[) are given respectively by Eq. (f3"3")l with G(p°) replaced by J(p°), and by 
Eqs. (|34H38|) with H(p°) replaced by J{p )- We immediately see that g\(x,pr) < fx(x,pr) and hi(x,px) < fx{x,Pr)- 
Since \P q \ < N q , this means that the 'trivial' inequalities \g\{x,pr)\ < fi{x,pr) and \h\(x,pr)\ < fi(x,pr) hold. 

Using the notation of the 'unintegrated' transverse moment of a TMD introduced in l|30p we obtain the following 
equalities among the 'bare' (flavourless) functions 

Ii(x,Pt) +9x(x,pr) = 2h 1 (x,p T ) (A2) 
Ix{x,Pt) - 9x(x,p T ) = -2hf^\x,p T ) (A3) 
gt£\x,prr = h$\x,pr? = ^ [fx{x, PT f ~ gx(x,p T ) 2 ^ (A4) 



iVc + 1 , r Nc-l 

—2— , ^ = — 2~ 

~~ ^ — ' Pd = ~a — ' and analog g± T , hx, h^ L , h^ T . (Al) 

6 o 
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The relations (|A2[ IA3[) among the bare distributions were discussed previously in various models [64], [6|| . It is 
important to notice that even if one assumes SU(6) symmetry, the relations (|A2[|A3|) among bare distributions do not 
need to imply relations among TMDs of definite flavour, though it is the case in the bag and constituent-quark models 
f65j] . But the spectator model of Ref. [63| provides a counter-example: there the bare TMDs satisfy (|A21 IA3|) . but 
the flavoured TMDs constructed from the do not. 

The equalities (|A2IIA4[) do not mean that the inequalities (l26l[29| are saturated. That would be the case, in SU(6), 
only for TMDs of s-quarks in A where N s = P s = 1, see [94| where the Soffcr bound 93| was discussed. For the 
nucleon in SU(6) we have \P q \ < N q , and the equalities (1A2HA4[) lead to real (never saturated) inequalities (j26l - [29|) . 
(We recall that T-odd distributions are absent in our approach.) 

Thus, we conlcude that the inequalities are manifestly satisfied in our approach — if one assumes SU(6) symmetry. 
If one does not, the positivity conditions (I26II29[) 'translate' into certain constraints among the covariant momentum 
distributions G(p°) and H(p°), see Ref. [84j. 



APPENDIX B: RELATIONS AMONG TMDS IN SU(6) 



In this Appendix we discuss relations among TMDs that are obtained in our approach under the assumptions of 
SU(6) symmetry and that at all TMDs are characterized in terms of the same covariant momentum distribution J(p°), 
see App.[5] From Eqs. (|A1| - |A4[) we obtain the following relations among TMDs with definite flavour 

%-ft(x,pT) + gKx,PT) = 2hl(x,p T ) (Bl) 

^-fI(x,PT)-9l(x,PT) = -2h^ 1)q (x, PT ) (B2) 

g^"( X ,p T f = h^ 1)q ( X ,p T r = ^^f«( X)PT f-g<l( X)PT )^ (B3) 

We remark that the relations (|Bli IB2|) hold in the bag (64| and constituent-quark [65l | model, but not in spectator 
models 63, 68] . Integrated versions of (|B1[) were discussed previously in [H, [94], HH| • 

The assumption of SU(6) symmetry by itself is a phenomenologically well-motivated concept especially in the 
valence- a; region, see [78| for a recent discussion in the context of TMDs. However, in our approach this is not yet 
a sufficient condition for the relations (|B1IIB3[) to be valid. In addition to SU(6) symmetry, we have to assume here 
that the covariant momentum distribution J(p°) appears in all TMDs. 

This fully supports the observation [ftjj . that SU(6) symmetry in a quark model alone is not a sufficient condition 
for this kind of relations to hold. Another SU(6) symmetric model, which in general does not support (|B1I IB3|) , is the 
spectator model of [63| — though upon an additional assumption (large-A^ limit) they hold there, too [64|. 



APPENDIX C: PROOF OF THE WW- TYPE RELATION EQ. (JUJ) 



In this Appendix we present two independent proofs of Eq. ([55)) . For the first proof, we use the notation of Eq. (22) 
in HI] to write the model expressions ([331 [37]) as 

^(^-^(xj + y^W + O^), hl{x) = ^ v _ 2{x) + (^) . (CI) 

Then, exploring the identity V^_ 1 (x) = f V_ 2 { x ) derived Eq. (24) of *86|, we obtain 

p^y _ *£, _ _ VLi(x) + * v , 2ix) +0 ^_ ^ (C2) 

which is equivalent to Eq. ([55)1 . 

For the second proof we show that hf^ q ( X ) and — x 2 f% h\(y) have the same Mellin moments. Notice that in the 
model all TMDs j q { X ) are well-behaving functions without singularities, have no support outside the region x £ [0, 1], 



14 



and have well-defined Mellin moments f dx x N j q (x) V N = 0, 1, 2, . . . Therefore Eq. (|55|) is equivalent to 

A^sjT'dx^^^+l^ftf^)) =o(j^) ■ (C3) 



Introducing the notation [dp 3 ] = ^-jy- H bf^ we write the expressions ([551 157)) for /if (a;) and hf^ q (x) as 



*!(*)= —]^\^-*)* 2 + °[m)- (C5) 

where p° — \p | (recall we neglect m). We insert the expressions (|C4[ IC5[) into Eq. (|C3[) . interchange the order of the 
integrations over x and p, introduce spherical coordinates such that p 1 — \p \ cos 9, and obtain 

An = / {d'p} \P\ N+2 ( - (1 - cos<+ 2 (1 + cosfl) + (1 + o(g) ■ (06) 



Now, our proof is completed because in (|C6[) the integral over z = cos 9 is 

iV + 3 / Li dz\ N + 3 
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